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Abstract 

We consider the constraints from proton decay and b-T unifica- 
tion in the minimal supersymmetric SU(5) grand unified theory with 
a 'visible' dynamical supersymmetry breaking sector. We show how 
the presence of vector-like messenger fields and the constrained super- 
particle mass spectrum affect the phenomenology of the model. We 
include the messenger fields in our renormalization group analysis be- 
tween the messenger scale (~ 100 TeV) and the GUT scale. We show 
that the simplest model of this type, a minimal SU(5) GUT with an 
additional 5+5 of messenger fields is excluded by the constraints from 
proton decay and b-T unification. 
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1 Introduction 



Supersymmetric grand unification is one of the viable possibilities for the 
physics that lies beyond the standard model 0. Interest in supersymmet- 
ric grand unified theories (SUSY GUTs) has been motivated primarily by 
two observations: (i) super symmetry eliminates the quadratic divergences in 
radiative corrections to the Higgs mass, which can destabilize the hierarchy 
between the GUT and electroweak scales 0, and (ii) gauge coupling unifi- 
cation can be achieved to considerable accuracy provided that superparticle 
masses are < 1 TeV jj. While this picture may turn out to be correct, it 
provides no explanation for why the supersymmetry breaking scale is so low 
(for example, in comparison to the Planck scale). In addition, it gives no 
explanation for the near degeneracy of the squark masses, necessary to avoid 
large flavor-changing neutral current (FCNC) effects. While supersymmetry 
can resolve (i) and (ii) above, the origin, scale and pattern of supersymme- 
try breaking masses required to produce a viable phenomenology remains a 
separate and important puzzle. 

Recently proposed models of dynamical supersymmetry breaking (DSB) 
provide possible solutions to these remaining problems ||, f| . If supersymme- 
try is broken by nonperturbative dynamics triggered when some asymptotically- 
free gauge coupling g(A) becomes large, then we would expect 

A^m exp — (1) 

bg 2 (m ) 

where b < is a beta function, and m is some high scale, like M Planck . 
The exponential suppression can account for a hierarchy between the SUSY- 
breaking scale and the Planck scale. Furthermore, the models that have been 
proposed so far employ a mechanism by which supersymmetry break- 
ing is transmitted to the ordinary particles through loop diagrams involving 
vector-like "messenger" fields, that carry electroweak quantum numbers, and 
ordinary gauge interactions. When the messenger fields feel SUSY breaking 
originating from the dynamical SUSY breaking sector of the theory, they 
transmit it to the ordinary squarks via these diagrams, which are flavor in- 
dependent ||, ||]. As as a result, the squarks of different generations remain 
degenerate, and the FCNC problem is naturally avoided. 

In this paper, we will comment on the GUT phenomenology of models 
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of this type. In the present context, GUT refers to the unification of the 
ordinary gauge groups of the standard model, but not to the unification 
of these groups with the additional gauge groups responsible for dynamical 
supersymmetry breaking. We restrict our discussion to the minimal SU(5) 
grand unified model || for definiteness. We will first argue that the existence 
of vector-like multiplets with electroweak quantum numbers at a scale ~ 
100 TeV is generic to any workable model. We will catalog the possible 
particle content of the messenger sector that is allowed by the requirement of 
perturbative gauge coupling unification, and we state the predicted squark 
and gaugino mass relations that follow in each case. Using this information, 
and including the messenger fields in our renormalization group analysis, we 
determine the bounds from the nonobservation of proton-decay, and from the 
requirement of b-r Yukawa unification. We conclude that a minimal SU(5) 
GUT with the simplest messenger sector possible is excluded by the lower 
bounds on the proton half-life. 



2 Messenger Sector 

The part of the messenger sector that is relevant to our analysis involves 
those fields which carry electroweak quantum numbers. These fields trans- 
mit SUSY breaking to the ordinary sector via loop diagrams involving elec- 
troweak gauge interactions. Since these diagrams have appeared in a number 
of places in the literature j|, |5|, [7| , we do not display them again here. In the 
model of Ref . , the relevant part of the messenger sector superpotential is 

W m = X D SDD + XiSll (2) 

where the fields have the quantum numbers D ~ (3, l)-i/3 D ~ (3, 1)1/3, 
/ ~ (l,2)_i/ 2 , and 7 ~ (1, 2)1/2 under the standard model gauge group, 
and S is a singlet chiral superfield. This particle content forms full SU(5) 
multiplets, 5 + 5, so that the apparent gauge unification at ~ 2 x 10 16 GeV 
is preserved. The vacuum expectation value (vev) of the scalar component 
of S, (S), determines the SU(3)xSU(2)xU(l) invariant masses in (0), while 
the vev of the F component of S, (F s ), parametrizes the degree of SUSY 
breaking. In Ref. || the 3-2 model |§ is assumed as the source of DSB, and 
the authors show that the remaining portion of the messenger superpotential 
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can be constructed so that DSB in the 3-2 sector generates vevs for both S 
and F s . 

We first would like to argue that the portion of the messenger superpo- 
tential given in (Q) is generic to a wide variety of realistic models in which 
SUSY breaking is transmitted to the ordinary sector via gauge interactions. 
If the messenger quarks and leptons are vector-like under the standard model 
gauge group, and we allow no dimensionful couplings in the superpotential, 
then the couplings in (0) will be present. The latter requirement is a philo- 
sophical one, namely, that all mass scales in the theory be generated via di- 
mensional transmutation. One might imagine constructing a model in which 
the messenger fields are chiral rather than vector-like under the standard 
model gauge group. The general problem of a model of this type is that 
radiatively-generated gaugino masses are too small. The one-loop diagram 
responsible for generating a gaugino mass necessarily involves chirality flips 
on the fermion and scalar lines. These chirality flips are proportional to the 
messenger fermion masses, rrif, and therefore the gaugino masses are of order 
(a/i^rrif/msusY- Since nif is of order the weak scale, the gauginos in this 
scenario are unacceptably light .[] 

One might also imagine that the messenger quarks and leptons are vector- 
like under the standard model gauge group, but carry nonstandard quantum 
numbers as well. In this case, it is likely that the additional gauge couplings 
would be perturbative. The only gauge group that is nonperturbative is the 
one in the DSB sector, and introducing additional particles that transform 
under it could lead to disaster in two ways: the vacuum structure of the 
theory may change so that supersymmetry is restored, or the multiplicity 
of messenger quarks and leptons may be too large to retain perturbative 
unification of the ordinary gauge couplings J3J. We know of no workable 
model in which the messenger quarks and leptons couple directly to a strongly 
interacting group from the DSB sector. If the messenger quarks and leptons 
couple to a nonstandard gauge group that is perturbative (like SU(3) in 
Ref. [|J), one might still worry that the conclusions of the two-loop proton 
decay analysis that we will present in the next section could be altered. While 
the standard model gauge couplings will run differently in this case, one can 

^We do not consider the possibility of light gluinos in this paper |J. 
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show that the quantities relevant to our analysis (e.g., the mismatch of gauge 
couplings at the GUT scale, mass ratios of particles at the messenger scale, 
etc.) will remain unaltered and our conclusions will remain the same. 

In what follows we assume minimal SU(5) unification, so that the gauge 
structure of the theory is SU(5) xG DSB . In Gdsb we include any nonstandard 
gauge groups that may be necessary for communicating SUSY breaking to the 
full messenger sector superpotential (e.g., the messenger hypercharge group 
discussed in ||). SU(5) gauge invariance implies that the messenger quarks 
and leptons form complete SU(5) representations. In the minimal case of a 
5+5 in the messenger sector, the messenger superpotential at the GUT scale 
has the form 

W m = XS55. (3) 

However, below the GUT scale, SU(5) is broken, and we recover the super- 
potential given in (|2]). The assumption of unification allows us to compute 
\e> and \i in terms of A, by running these couplings down to the messenger 
scale. Once the messenger scale has been specified, threshold corrections at 
this scale are calculable, and can be included in our renormalization group 
analysis without introducing any additional uncertainty. This is true for 
representations larger than 5+5 as well. 

Next, we must specify what other SU(5) representations are allowed in the 
messenger sector. Introducing additional SU(5) multiplets preserves gauge 
unification, but the gauge coupling at the GUT scale ot^imcuT) increases 
as we add additional multiplets. If we require that a^rriGUT) remains per- 
turbative, then we may add 1, 2, 3 or 4 (5+5) pairs, or a single (10+10) 
pair, or (5+5) + (10+10) to the particle content of the minimal SU(5) GUT. 
Additional 5s or 10s, or larger SU(5) representations will render a^mGUT) 
nonperturbative |L0|. 

The messenger sector fields not only affect the renormalization group 
analysis between the messenger scale A and the GUT scale (including the 
calculable threshold corrections at the scale A) but also constrain the thresh- 
old corrections at the weak scale. As presented in Ref. |5]], the radiatively- 
generated gaugino and squark masses, rrii and rh, assuming a 5+5 in the 
messenger sector, are given by 

m- = (4) 

1 167T2 (S) ' 1 ] 
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i6vr 2 ; 

where C F is 3/4 for SU(2) doublets, 4/3 for SU(3) triplets, and (3/5)F 2 for 
U(l). The quantity (Fs)/(S)= mo has dimensions of mass, and is of the 
same order as the messenger scale A. In a specific model, mo is calculable 
by minimizing the messenger sector potential. Once is fixed, the parti- 
cle content of the messenger sector completely determines the gaugino and 
squark masses at the scale A; these can subsequently be run down to the 
weak scale. We explain how we fix the precise value of m in the following 
section. If the messenger sector consists of n 5 5+5 pairs, then ([§) and @ 
both scale as 715, while the ratio rrii/m scales as y/n*,. In the case of a 10+10 
pair, we obtain the same result for the gaugino and squark masses as the 
case where n 5 = 3. These observation will be useful in our discussion of the 
proton decay bound in the following section. 



3 Proton Decay Analysis 

Our analysis of the proton decay constraints is similar in spirit to that of 
Refs. [p+| , |T2|| . By including threshold corrections at the GUT scale, we can 
determine the largest color-triplet Higgs mass that is consistent with gauge 
coupling unification. We can then constrain the remaining free parameters 
involved in the proton decay matrix element using the current lower bounds 
on the proton half-life. Our algorithm is as follows: 

i. We fix 115 or uio, the number of 5+5 or 10+10 pairs in the messenger 
sector. The messenger sector scale is set at 100 TeV, to insure that the 
superparticle masses are of order the weak scale (recall equations (f3j) and 
above) . 

%%. We define the GUT scale M GUT as the scale where the SU(2) and 
U(l) gauge couplings unify. We use the input values af 1 (m^) = 58.96 + 0.05 
and aj" 1 ^) = 29.63 ± 0.05, that follow from Ref. [|13|]. We first run them 
up to m top = 176 GeV using the standard model renormalization group 
equations (RGEs). We then numerically solve the two- loop RGEs with the 
super symmetric particle content, taking into account the change in the one- 
and two-loop beta functions as we cross the messenger scale. The RGEs and 
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beta functions are provided in the appendix. 

Hi. Given an input value of (mz), we determine the mismatch between 
«3 1 (Mgut), and (Mqut), the GUT coupling determined in step ii. We 
ascribe this mismatch to the sum effects of threshold corrections at the weak 
scale, the messenger scale, and the GUT scale. Using the one-loop expressions 



for threshold corrections given in Ref . [llj , we find that the mismatch Aa 3 
a^ 1 (M GUT ) - a^ 1 (M GUT ) is given by 

1 



47T 



5 M G u T 



m 



h 



+ 4 In In 

TTlyj 5 fTlfop 



12 m D 
+ — n 5 In 

5 rriL 



(6) 



The first term gives the threshold correction at the GUT scale as a function 



of the color-triplet Higgs mass M# c . The next two terms give the largest [11 



threshold corrections at the weak scale, depending on the gluino, wino, and 
higgsino masses. We checked the threshold corrections from scalars can be 
safely neglected in this analysis. The third term gives the threshold correction 
at the messenger scale, from the splitting of the original n 5 5+5 pairs. In the 
case of a 10+10 pair, we have a different particle content at the messenger 
scale, and we should make the replacements 



12 m D 
— n 5 log 

5 wil 



6 mg 18 mg 
--log -log 

5 niE 5 mu 



(7) 



where the new fields have the quantum number assignments Q ~ (3, 2)i/ 6 , 
andtf~(3,l)_ 1/3 . 
Note that most of the variables in (0) and (|7|) can be estimated reli- 
ably enough, making it possible to place an upper bound on Mh c - We take 
fng/n^w — ~ 3.5, and set m~ h = 1 TeV to maximize M# c . The ratio 

mg/m^ can be computed by running the Yukawa couplings in (0) between 
the GUT scale and the messenger scale. We use a one-loop estimate of this 
ratio which takes into account the effect of the gauge interactions on the 
running. We find 



m D /m L 



a 5 1 


1 1 

6 TT+"5 


a 5 1 


3 1 

2 l+ri5 


a 5 1 


8 1 
3 -3+n 5 






|V(A)_ 









(8) 



for 725 < 3. For = 3 (the case where a.3 no longer runs at the one- loop 
level above the messenger scale), one makes the substitution 



a 5 1 


8 1 

3 — 3+nr, 


f A " 


4 c«3( A ) 

3 7T 


[a 3 - 1 (A)_ 




-M GUT . 





(9) 
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in eq. Finally, the ratios that are relevant in the 10+10 case are given 
by 



m Q /m E 



"5 1 


288 


a 5 x 
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\ A 1 






k _1 (A)_ 




-M GUT . 



4 t»3 (A) 

3 TV 



mg/mu 





a 5 1 


15 
288 


a 5 1 


3 
8 








|V(A)_ 





(10) 

(11) 



For any input value of 03(1712), we determine the mismatch Aaj$ , and 
the maximum Mh c that follows from (^j) or (^). Note that in determining this 
upper bound, we take into account the uncertainties in Mq UT and a 5 (M GUT ) 
that follow from the experimental uncertainties in ot\{mz) and « 2 (m2) at 
90% confidence level. The results are shown in Figure 1. In the case where 
we ignore threshold corrections at the messenger scale, Figure la, we see that 
the addition of 5+5 pairs tends to weaken the bound on M# c ; in the 10+10 
case the bound is strengthened. However, we see that in the final result, 
Figure lb, the upper bound on Mjj c is not much different from the minimal 
case. This is due to an accidental cancellation between two different effects, 
namely the two-loop contribution to the gauge coupling evolution due to the 
additional fields, and the threshold correction at the messenger scale. 

iv. Using our results shown in Figure 1, we can study the bounds from 
the proton life-time. Consider the mode n 
given in Ref . [ 1 1 ] as 



K u^. The nucleon life-time is 



Tin 



K V^) = 3.9 x 10 di yrs. x 



0.003GeV 3 0.67 sin2/5 M Hc 



TeV 



-1 



C A s l + ^KFGeV/(M) + /Ke) " (12) 
Here ( parametrizes the uncertainty in the hadronic matrix element, and is 
estimated to be between 0.003 GeV 3 to 0.03 GeV 3 0. We set ( = 0.003 to 
be the most conservative. As represents the short distance renormalization 
of the proton-decay operators, i.e. from running the Yukawa couplings up 
to the GUT scale, and then running the dimension-5 operators generated 
by integrating out the color-triplet Higgs, back down to the weak scale. We 
have checked that the effect of the messenger sector fields on the numerical 
value of is negligible, and set As ~ 0.67. Note that this is an extremely 
conservative assumption. The effect of going to = 2, for example, reduces 
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As to 0.63. However, we have not included the effect of the top quark Yukawa 
coupling in the running which can enhance As (and strengthen the resulting 



proton decay bound) by as much as a factor of ~ 3 []TT| . We simply set 
As = 0.67 in our numerical analysis as a conservative estimate. / is the 
"triangle" function obtained by dressing the dimension-5 supersymmetric 
operators with winos to produce dimension-6 four-fermion operators that 



are responsible for the decay |Lo], lib] . It is given by 

/(«, d) = -^— 2 ( In 4 - -^hf ^ 4 V ( 13) 

The squark, slepton and wino masses fed into the triangle function are those 
given by the generalization of ([|) and (|) to the case of n 5 5+5 pairs or 
nio 10+10 pairs, as we discussed earlier. The precise values that we choose 
depend on the value of the parameter m , which is of order the messenger 
scale. We choose mo such that the squark masses are normalized 1 TeV 
at the weak scale, and the gaugino masses are fixed by (f§). We make this 
choice because the proton-decay bounds are weakest for the heaviest squark 
masses. However, we do not allow squark masses greater than 1 TeV in the 
interest of naturalness. Finally, the quantity 1 + y tK parametrizes possible 
interference between the diagrams involving c and t exchange |T6| |. While 
the bound from the Kv mode can be made arbitrarily weak by assuming 
destructive interference between diagrams (i.e. y = —1.0), one cannot 



simultaneously weaken the bounds from the other possible decay modes [III 
For example, for |1 + y t x\ < 0.4, the mode n — > ti°V^ becomes dominant, 
giving a comparable decay rate. Thus, we use the Kv mode with |1 + ytx\ = 
0.4 to obtain a characteristic bound on the parameter space. 

With all the parameters fixed as described above, we obtain a bound on 
the quantity sin 2/3, or alternately tan/?, the ratio of up- to down-type Higgs 
vevs, for each value of a^mz) that we input in step ii. The results are shown 
in Figures 2a through 2c. The area of the «3-tan (3 plane that is consistent 
with the proton decay constraints is the region above the solid line labeled 
proton decay. The region between the horizontal dashed lines show the area 
allowed by the two standard deviation uncertainty in the experimental value 
of a^imz) measured at LEP 0.116 ± 0.005 fl3| . It is clear that in going 
from n 5 = n 10 = to n 5 = 1, these regions no longer overlap. In the 
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^5 = ^io = case, we have chosen the most conservative set of parameters 
possible, namely, we have set all the squark and slepton masses fa — 1 TeV, 
and the wino mass = 45 GeV. However, in the 715 = 1 case, with squark 
masses normalized to 1 TeV, the slepton and gaugino masses are predicted 
from (f|) and fl5|). Because we have greater predictivity in this case, we 
can no longer choose the most favorable parameter set, and the bound is 
strengthened. This effect is large enough to overwhelm the competing effect 
of the slightly weaker bound on the color triplet Higgs mass that we found 
in Figure 1. As we pointed out earlier, the ratio of gaugino to squark mass 
scales at -y/fis, so as we go to larger values of n 5 , we obtain an even less 
favorable set of parameters. Our proton decay bound stays more or less the 
same however, due to the competition between this effect and the weaker 
bound on the color triplet Higgs mass. Finally, we obtain a slightly tighter 
bound in the 10+10 case, Figure 2c, as a result of the slightly tighter upper 
bound on Mh c shown in Figure 1. For each case in which a messenger sector 
is present, the proton decay region never intersects with the allowed region 
for a 3 (m z ). 

4 b-r unification 

For completeness, we have also determined the region of the ct3-tan /3 plane 
that is consistent with b-r unification. Here we work only to the one-loop 
level. In our proton-decay analysis, the bound on Mh c followed from thresh- 
old corrections, so we needed to include all other effects of equal importance. 
This necessitated the two-loop analysis. In the case of b-r unification, how- 
ever, we do not need the higher level of accuracy as we will explain below, 
so we worked at one loop. 

Our algorithm is straightforward. For a given choice of 03 and tan/3, we 
determine the top, bottom, and tau Yukawa couplings. We run these Yukawa 
couplings up to the scale m top using the standard model renormalization 
group equations. Above m top we run these Yukawa couplings using the one- 
loop supersymmetric renormalization group equations up to the GUT scale 
[ p| ] . We impose two conditions to determine whether we have a valid solution 
with b-r unification: (a) we require A;, and A T to be within 0.5% of each other 
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at the GUT scale (b) we require A&, A T and Aj op to be less than 2. Condition 
(b) is imposed so that the couplings do not blow up below 10 times the GUT 
scale. If the cutoff scale is lower than this, higer- dimension operators can 
generate corrections to the Yukawa unification at more than the 10% level 
[ |P?jj . The messenger sector particles alter the analysis through their effect 
on the running of the gauge couplings, which in turn enter into the one-loop 
RGEs for the Yukawa couplings. These renormalization group equations are 
provided in the appendix. The requirement of unification to within 0.5% at 
the GUT scale is somewhat arbitrary. Threshold corrections and the effects of 
^Planck supressed operators could in principle account for a larger mismatch 
between the Yukawa couplings at the GUT scale. Our results should therefore 
be considered qualitative. Unlike our proton decay analysis, we do not do 
a two-loop analysis including threshold corrections at each of the relevant 
scales. The presence of M P i anck supressed operators can give important GUT 
scale threshold corrections in the small tan/? region [17], rendering the extra 
accuracy of such an analysis meaningless. 

Our results are shown in Figures 2a through 2c. The crescent shaped 
region is generated by allowing the 6-quark MS mass to vary between 4.1 
and 4.5 GeV, the range suggested by the QCD sum rule analysis Note 
that as we increase n*, or n\o, the b-r region moves towards smaller values 
of «3. While this is not inconsistent with the measured value of «3, the 
distance between the b-r region and the area allowed by the proton decay 
bounds increases monotonically with or n w .f\ 



5 Conclusions 

We have shown that the minimal supersymmetric SU(5) model cannot be em- 
bedded successfully within the simplest type of scenario suggested in Ref . [|] . 

■i-Note that the low tan/3 'cusp' of the crescent would extend into the proton decay 
region for allowed as{mz) in the minimal case (71.5 = nio = 0) had we imposed the less 
stringent constraint X top < 3.3 Q on the acceptable size of the Yukawa couplings at the 
GUT scale. Even in this case, however, the three-way overlap between proton decay, b-r 
unification and a^{mz), does not persist when the messenger fields are present. If we had 
chosen a smaller value of mt op = 150 GeV, on the other hand, the b-r region would extend 
to smaller tan/3, but not to larger a^,{mz), and our conclusions would remain the same. 
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The effect of messenger sector particles on the renormalization group analy- 
sis, as well as the predicted gaugino and squark mass ratios that follow from 
the messenger sector particle content lead to a conflict with the lower bound 
on the proton lifetime. In addition, the region of parameter space preferred 
for b-T unification moves farther away from the region preferred by the pro- 
ton decay bounds as the size of the messenger sector is increased. We must 
emphasize that this is only a mild obstacle to the type of scenario proposed 
in Ref. ||. Non-minimal GUTs can be constructed which evade the proton 
decay bound p0| . It is nonetheless interesting that we can exclude what is 



perhaps the simplest models of this type. 
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A Appendix 

The two-loop supersymmetric renormalization group equation for the gauge 
couplings that we use in our proton decay analysis is 

"t = is? i « f + (is?)g^ (14) 

where the beta functions are given by 
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Here n g , n^, n 5 and n w are the number of generations (3), higgs doublets 
(2), messenger sector 5+5 and 10+10 pairs, respectively. 

The one-loop renormalization group equations for the top, bottom, and 
tau Yukawa couplings used in our analysis of b-r Yukawa unification are given 
by @ 
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Figure Captions 



Fig. 1 (a) Upper bound on the color-triplet Higgs mass M Hc as a func- 
tion of a.z{mz), ignoring threshold corrections (i.e. mass splittings) at the 
messenger scale. The solid line is the minimal SU(5) result, the dashed line 
is the case of three 5+5 pairs in the messenger sector, and the dotted line 
is the case of one 10+10 pair, (b) The complete result, including threshold 
corrections at the messenger scale. 

Fig. 2(a) Preferred regions of the tan(f3)-a 3 (m z ) plane, in minimal 
SU(5) unification, for m top = 176 GeV. The region above the proton-decay 
line is allowed by the lower bound on the proton lifetime (see the text) at 
90% confidence level, while the region within the crescent shaped contour is 
preferred by the constraint of b-r Yukawa unification with mbira^) = 4.1- 
4.5 GeV. The horizontal band shows the experimentally allowed range of 
OLz{mz) = 0.116 ± 0.005 at two standard deviations, (b) Same as (a) for 
n 5 — 1, (c) Same as (a) for n 5 = 3. The dashed line shows the result for 
n w = 1 and n 5 = 0; the b-r region remains the same for = 3 or n w = 1. 
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